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We  evaluate  the  theoretical  performance  of  a  point-source  integrating-cavity  absorption  meter  (PSICAM) 
with  Monte  Carlo  simulations  and  a  sensitivity  analysis.  We  quantify  the  scattering  errors,  verifying 
that  they  are  negligible  for  most  ocean  optics  applications.  Although  the  PSICAM  detector  response  is 
highly  sensitive  to  the  value  of  the  wall  reflectivity,  the  absorption  of  an  unknown  fluid  can  be  accurately 
determined  with  a  PSICAM  if  appropriate  reference  solution(s)  are  chosen.  We  also  quantify  the  error 
that  results  if  the  source  is  not  perfectly  isotropic,  finding  that  moderate  amounts  of  source  anisotropy  can 
be  tolerated  provided  that  the  detector  is  properly  located  with  respect  to  the  source. 
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1.  Introduction 

The  measurement  of  light  absorption  by  seawater 
and  its  components  is  of  great  interest  for  oceanog¬ 
raphy  and  for  optical  remote  sensing  of  the  ocean. 
Biological  oceanographers  measure  the  absorption  of 
light  by  phytoplankton  to  quantify  photosynthesis; 
chemical  oceanographers  are  concerned  with  the  ab¬ 
sorption  of  light  by  organic  compounds  and  the  photo¬ 
oxidation  of  those  compounds;  and  physical 
oceanographers  measure  the  absorption  of  energy  by 
surface  waters  and  its  role  in  the  heating  of  the  sur¬ 
face  layer  and  the  formation  of  the  thermocline.  Op¬ 
tical  oceanographers  are  concerned  about  all  these 
processes  as  well  as  use  of  spectral  absorption  mea¬ 
surements  to  identify  and  quantify  optically  active 
components  in  the  water. 

There  are  several  methods  for  determining  the  ab¬ 
sorption  coefficient  of  ocean  waters.  The  most 
widely  used  approach  is  the  quantitative  filter  tech¬ 
nique,1  in  which  water  samples  are  concentrated  onto 
glass-fiber  filters  and  the  absorption  of  the  resulting 
filters  and  filtrate  are  measured  separately  in  a  spec¬ 
trophotometer.  This  method  requires  corrections 
for  scattering  errors  associated  with  the  filter  and  for 
scattering  errors  in  the  spectrophotometric  measure¬ 
ment  of  the  filtrate.  Because  of  its  convenience  com¬ 
pared  with  the  quantitative  filter  technique,  use  of  in 
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situ  absorption  meters,  in  particular  the  ac-9  (WET 
Labs,  Inc.,  Philomath,  Oregon),2-3  is  becoming  in¬ 
creasingly  popular.  These  instruments  measure  the 
attenuation  of  an  active  light  source  through  a  vol¬ 
ume  of  seawater,  but,  like  a  spectrophotometer,  must 
be  corrected  for  scattering  errors.  A  third  approach 
to  determining  the  absorption  coefficient  of  seawater 
is  to  measure  the  natural  light  within  the  water  col¬ 
umn  and  to  invert  the  radiative  transfer  equation  to 
calculate  the  values  of  the  water’s  optical  properties 
that  are  consistent  with  the  observed  light  measure¬ 
ments.4-5  The  primary  difficulty  with  this  approach 
is  the  need  to  remove  numerically  or  analytically  the 
effects  of  the  boundary  conditions  on  the  light  field. 

An  alternative  to  these  three  approaches  is  use  of 
an  integrating-cavity  absorption  meter  (ICAM), 
which  is  the  focus  of  this  research.  The  motivation 
of  this  approach  is  to  obtain  a  measurement  that  is 
insensitive  to  scattering.  The  water  sample  is  ac¬ 
tively  illuminated  in  the  cavity  of  the  ICAM,  the  walls 
of  which  are  made  of  a  highly  reflective  material  so 
that  the  resulting  light  field  within  the  cavity  is 
nearly  uniform.  The  absorption  of  the  sample  is 
taken  to  be  the  power  input  to  the  ICAM  minus  either 
the  measured  or  the  calculated  losses  at  the  cavity 
wall. 

Elterman6  suggested  the  use  of  a  highly  reflective 
Lambertian  cavity  to  obtain  a  scattering-insensitive 
measurement  of  the  absorption  coefficient  of  a  solid 
material.  A  small  sample  of  the  material  was  placed 
in  an  otherwise  empty  cavity,  and  the  light  source 
and  detector  were  placed  at  the  cavity  walls.  Fry 
and  Kattawar7-8  later  adapted  this  approach  for  de¬ 
termining  the  absorption  coefficient  of  seawater. 
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Their  configuration  consists  of  a  cavity  completely 
filled  with  the  water  sample  and  a  second  integrating 
cavity  surrounding  the  first.  An  isotropic  light  field 
is  generated  between  the  two  cavities  that  diffuses 
into  the  inner  cavity,  providing  to  the  inner  cavity  a 
diffuse  light  source  that  is  uniform  over  the  cavity 
walls.  Kirk9  proposed  that  the  two  cavities  be  made 
spherical  and  concentric,  making  it  possible  to  model 
the  relationship  analytically  between  the  absorption 
coefficient  and  the  ICAM  detector  response. 

Kirk10  later  proposed  a  different  arrangement  in 
which  the  spherical  cavity  containing  the  water  sam¬ 
ple  is  illuminated  with  an  isotropic  source  at  the 
center  of  the  sphere.  The  primary  advantage  of  this 
point-source  integrating-cavity  absorption  meter 
(PSICAM)  over  the  concentric  cavity  design  is  the 
elimination  of  the  outer  cavity,  which  is  difficult  to 
construct.  Kirk10  provides  equations  for  obtaining 
the  absorption  coefficient  with  a  PSICAM  and  in¬ 
cludes  some  analysis  of  the  design. 

Here  we  extend  Kirk’s  PSICAM  analysis  in  prepa¬ 
ration  for  the  construction  of  a  prototype  instrument 
that  will  be  used  at  the  U.S.  Naval  Research  Labora¬ 
tory  for  the  investigation  of  coastal  ocean  waters. 
Kirk’s9-10  equations  for  the  operation  of  a  PSICAM, 
which  ignore  scattering  effects,  are  reviewed  in  Section 
2.  The  errors  in  these  equations  that  are  due  to  scat¬ 
tering  are  quantified  in  Section  3,  and  a  sensitivity 
analysis  of  the  equations  is  provided  in  Section  4.  In 
Section  5  we  investigate  the  error  that  is  introduced 
when  the  light  source  is  not  perfectly  isotropic.  A 
summary  of  our  conclusions  is  provided  in  Section  6. 


Therefore 


P0A(  1  -  P PsB) 

p0B(l-ppsA)- 


(3) 


For  a  central  point  source  the  value  of  P0  for  a  non¬ 
scattering  solution  is10 


P0(a,  r )  =  exp(-ar),  (4) 

where  a  is  the  absorption  coefficient  of  the  sample 
and  r  is  the  inner  radius  of  the  PSICAM  cavity. 
However,  if  the  source  is  generated  with  a  diffusing 
sphere  that  has  a  radius  rs  that  is  significantly  large 
compared  to  the  cavity  radius,  then  we  should  take  P0 
to  be 


P0{a,  r0)  =  exp(-ar0),  (5) 

where  r0  =  r  —  rs.  For  diffuse  light  leaving  the 
spherical  cavity  wall  the  value  of  Ps  for  a  nonscatter¬ 
ing  fluid  is9 


Ps(a,  r )  =  — [1  -  exp(— 2ar)(2ar  +  1)].  (6) 

2a  r 

The  size  of  the  source  is  unlikely  to  be  significant  for 
Ps  provided  that  the  source  material  is  nonabsorbing. 
For  very  small  values  of  ( ar ),  for  which  Ps  »  1,  Eq.  (6) 
is  difficult  to  evaluate;  however,  the  power  series9 

Ps(a,  r)  =  1  -  -  (ar)  +  (ar)2  -  (ar)3 

3  15 


2.  Basic  Equations 

The  outwardly  directed  irradiance  F0  at  the  inner 
wall  of  a  perfectly  symmetrical  PSICAM  is  propor¬ 
tional  to  the  average  number  of  times  a  photon  col¬ 
lides  with  the  wall,  Nc,  before  being  absorbed  either 
by  the  fluid  in  the  cavity  or  by  the  cavity  wall.9-10 
The  transmittance  Tab  that  is  the  ratio  of  the  mea¬ 
sured  values  of  F0  when  the  cavity  is  filled  with  sam¬ 
ples  A  and  B,  respectively,  is  therefore 


Tar  — 


Fl 

F0B 


N l 
Ncb 


(1) 


Note  that  the  effects  of  the  source  strength  and  de¬ 
tector  gain  on  F0  cancel  out  when  forming  the  ratio 
Tab ■  The  value  of  Nc  depends  on  three  quantities: 
(1)  the  probability  P0  that  a  photon  leaving  the  source 
at  the  center  reaches  the  wall  (i.e.,  without  being 
absorbed  by  the  fluid),  (2)  the  wall  reflectivity  p,  and 
(3)  the  probability  Ps  that  a  photon  leaving  the  wall 
will  return  to  the  wall.  The  value  of  Nc  equals  the 
number  of  photons  reaching  the  wall  directly  from 
the  source  plus  the  number  of  photons  colliding  with 
the  wall  for  a  second  time  plus  the  number  colliding 
with  the  wall  for  a  third  time,  etc.: 


+  •  •  •  2(-2ar)" 
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converges  rapidly  for  very  small  values  of  (ar). 
From  Eqs.  (3),  (5),  and  (6),  the  transmittance  of  a 
nonscattering  fluid  in  a  PSICAM  with  a  diffuse  cavity 
wall  is 


Tab  — 


exp[— r 0(aA  -  aB)] 


1  -  p Ps(aB,  r ) 
1  -  p Ps(aA,  r) 


(8) 


with  Ps(a,  r)  given  by  Eq.  (6). 

Kirk10  indicates  that  both  P0  and  Ps  are  insensitive 
to  the  presence  of  scattering  for  most  practical  ocean 
optics  purposes  and  that  Eqs.  (6)  and  (8)  can  therefore 
be  applied  to  seawater  samples  even  though  scatter¬ 
ing  was  ignored  in  their  derivation.  In  Section  3  we 
quantify  the  magnitude  of  the  scattering  error  in 
these  equations  for  a  large  range  of  physical  param¬ 
eters. 

Equation  (8)  can  be  used  in  two  ways.  First,  when 
the  absorption  coefficients  of  samples  A  and  B  are 
known  the  wall  reflectivity  can  be  determined  from 
the  measurement  of  with  the  solution10  of  Eq.  (8): 


NC  =  P0  +  PoPPs  +  PoP2Ps2  +  .  .  .  =  P0  2  (P PsT 

n- 0 

=  P0/(l  -  p Ps).  (2) 


_ _ Tab  exp(-aBr0)  -  exp(-qAr0) _ 

P  Tab  exp (-aBr0)Ps(aA,  r)  -  exp (-aAr0)Ps(aB,  r) ' 

(9) 
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Second,  from  Eq.  (8)  we  can  express  the  transmit¬ 
tance  T  of  a  fluid  of  interest  with  unknown  absorption 
coefficient  a  with  respect  to  a  reference  sample  R  with 
known  absorption  coefficient  aref: 


T  =  exp[— r 0(a  -  aref)] 


1  -  pPs(aref,  r) 
1  -  p Ps(a,  r ) 


(10) 


The  value  of  a  can  therefore  be  determined  from  the 
measurement  of  T  with  an  iterative  solution  of  Eq. 
(10). 

In  summary,  we  can  determine  the  unknown  ab¬ 
sorption  coefficient  a  of  a  sample  of  interest  by  (1) 
measuring  the  irradiance  values  at  the  PSICAM  cav¬ 
ity  wall  with  the  cavity  filled  with  solutions  A  and  B 
(of  known  absorption  coefficients  aA  and  aB)\  (2)  form¬ 
ing  the  ratio  ^AB  from  these  two  measurements  [Eq. 
(1)];  (3)  computing  p  with  Eq.  (9);  (4)  measuring  the 
irradiance  at  the  PSICAM  cavity  wall  both  with  the 
cavity  filled  with  the  sample  of  interest  (with  un¬ 
known  absorption  coefficient  a)  and  with  the  cavity 
filled  with  reference  solution  R  (with  known  absorp¬ 
tion  coefficient  aref);  (5)  forming  the  ratio  T  of  these 
two  measurements;  and  (6)  solving  Eq.  (10)  for  a. 
Note  that  the  required  number  of  reference  solutions 
with  the  known  absorption  coefficient  can  be  reduced 
from  three  (R,  A,  and  B)  to  two  (R  and  A)  if  solutions 
B  and  R  are  taken  to  be  the  same.  Furthermore,  if 
the  value  of  p  is  known  from  direct  measurement, 
only  steps  (4)-(6)  are  needed  and  we  require  only  one 
reference  solution  (R). 


3.  Scattering  Effects 

Equations  (5)  and  (6)  for  the  probabilities  of  photon 
survival  P0  and  Ps  were  derived  for  nonscattering 
absorbers.  In  Monte  Carlo  simulations,  Kirk9-10 
found  P0  and  Ps  for  seawater  to  be  insensitive  to  the 
value  of  b.  However,  in  general,  the  presence  of 
scattering  decreases  the  value  of  P0  by  increasing  the 
effective  path  length  of  each  photon  from  the  source 
to  the  wall.  Scattering  also  alters  the  effective  path 
length  of  a  photon  as  it  travels  from  the  wall  back  to 
the  wall.  Therefore  there  must  be  limits  on  the  val¬ 
ues  of  a,  b,  and  r  for  which  the  application  of  Eqs.  (5) 
and  (6)  is  valid. 

We  examined  the  effect  of  scattering  on  the  values 
of  P0,  Ps,  and  Nc  with  Monte  Carlo  simulations. 
These  simulations  were  performed  for  both  isotropic 
scattering  and  for  scattering  given  by  the  average  of 
Petzold’s  seawater  measurements  in  San  Diego  Har¬ 
bor.11  The  Petzold  scattering  phase  function  is 
highly  peaked  in  the  forward  direction,  with  back- 
scattering  representing  only  approximately  2%  of  the 
total  scattering  (as  opposed  to  50%  for  isotropic  scat¬ 
tering).  The  values  of  P0  and  Ps  were  computed  sep¬ 
arately,  and  Eq.  (2)  was  used  to  determine  Nc  for 
specific  values  of  p.  The  key  equations  for  the  Monte 
Carlo  simulations  are  given  in  Appendix  A. 

Because  the  value  of  P0  decreases  with  increasing 
values  of  b,  Eq.  (5)  overestimates  the  true  value  of  P0. 
The  percentage  error  in  Eq.  (5)  was  found  to  increase 
approximately  linearly  with  respect  to  b  and  to  be 


1  10  100 
b  mi l) 


Fig.  1.  Percentage  scattering  errors  in  P0  and  1/(1  -  p Ps)  for 
isotropic  and  Petzold  scattering  when  the  probabilities  of  source- 
to-wall  photon  survival  (P0)  and  wall-to-wall  photon  survival  (Ps) 
are  computed  with  Eqs.  (5)  and  (6)  with  a  =  1.0  m_1  and  r  =  r0  = 
0.05  m. 


greatest  for  large  values  of  a  and  r0.  The  error  for 
Petzold  scattering  was  much  smaller  than  that  for 
isotropic  scattering  because  for  Petzold  scattering  the 
photons  are  scattered  predominantly  in  the  near¬ 
forward  direction.  For  example,  the  percentage  er¬ 
ror  in  P0  determined  from  Eq.  (5)  is  shown  versus  b 
(log  scale)  in  Fig.  1  for  radius  r0  =  r  =  0.05  m,  wall 
reflectivity  p  =  0.99,  and  absorption  coefficient  a  =  1 
m  '.  For  isotropic  scattering,  the  error  in  Eq.  (5)  is 
approximately  0.9%  and  11%  for  b  =  10  m  1  and  b  = 
100  m  ',  respectively,  whereas  for  Petzold  scattering 
the  error  in  Eq.  (5)  is  only  approximately  0.04%  and 
0.6%,  respectively.  Fora  =  0.1  m  ',  the  errors  were 
only  approximately  one  tenth  of  those  for  a  =  1  m  1 . 

When  a  photon  traveling  from  the  wall  back  to  the 
wall  is  scattered,  its  effective  path  length  may  be 
smaller  or  larger  than  it  would  have  been  if  unscat¬ 
tered,  depending  on  where,  and  into  which  direction, 
it  is  scattered.  For  small  values  of  a,  b,  and  r  it  was 
found  that  the  overall  (statistical)  scattering  effect  is 
negligible.  For  large  values  of  a,  b,  or  r,  however,  it 
is  more  likely  that  a  photon  will  be  scattered  back  to 
the  wall  near  where  it  left,  thus  increasing  the  value 
of  Ps,  than  it  is  that  its  path  length  will  be  increased 
by  scattering.  Equation  (6)  therefore  underesti¬ 
mates  the  true  value  of  Ps  for  large  values  of  a,  b,  or 
r.  Even  if  the  scattering  error  in  Ps  appears  to  be 
insignificantly  small,  it  may  nonetheless  have  a  sig¬ 
nificant  effect  on  the  total  PSICAM  response  through 
the  term  1/(1  —  p Ps).  The  percentage  error  in 
1/(1  —  pPs)  is  shown  in  Fig.  1  for  r  =  0.05  m,  p  =  0.99, 
and  a  =  1  m  '.  Its  magnitude  increases  approxi¬ 
mately  linearly  with  increasing  values  of  b  but  is 
approximately  one  third  that  of  the  error  in  P0. 

Because  scattering  affects  P0  and  Ps  in  opposite 
ways,  the  percentage  error  in  F0  (which  equals  the 
percentage  error  in  Nc)  is  smaller  than  that  in  either 
P0  or  Ps.  The  error  in  the  PSICAM  response  pre¬ 
dicted  by  Eqs.  (2),  (4),  and  (6)  is  shown  in  Fig.  2  for 
isotropic  scattering  and  in  Fig.  3  for  Petzold  scatter¬ 
ing.  The  scattering  effect  on  F0  was  found  to  be 
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Fig.  2.  Percentage  error  in  the  PSICAM  response  predicted  by 
Eqs.  (2),  (4),  and  (6)  that  is  due  to  the  presence  of  isotropic  scat¬ 
tering  for  r  =  0.05  m  and  p  =  0.99. 


insignificant  for  most  oceanic  problems;  even  for  the 
extremely  large  value  of  b  =  100,  the  error  in  F0  for 
Petzold  scattering  is  less  than  0.35%  when  a  <  lm-1. 
However,  the  scattering  effects  may  be  important  for 
the  infrared  wavelengths  (for  which  a  >  1  m  ')  or  for 
other  PSICAM  applications  that  have  a  less  anisotro¬ 
pic  scattering  phase  function  than  water  or  have  sig¬ 
nificantly  larger  values  of  a,  r,  or  b  than  those  used  in 
Fig- 3. 

4.  Sensitivity  Analysis 

The  determination  of  the  absorption  coefficient  of  a 
fluid  sample  with  Eq.  (10)  requires  several  measure¬ 
ments,  the  uncertainties  of  which  result  in  a  total 
uncertainty  in  the  estimate  of  a.  Here  we  examine 
the  sensitivities  of  the  determined  value  of  a  to  these 
measurements.  We  first  conduct  this  analysis  for 
the  mode  of  PSICAM  operation  in  which  the  cavity- 
wall  reflectivity  p  is  assumed  to  be  known  from  direct 
measurement.  We  then  focus  on  the  mode  of  oper¬ 
ation  in  which  p  is  instead  determined  from  PSICAM 
measurements  of  two  nonscattering  samples  of 
known  absorption  coefficient  values  [i.e.,  with  Eq. 

(9)]- 


Fig.  3.  Percentage  error  in  the  PSICAM  response  predicted  by 
Eqs.  (2),  (4),  and  (6)  that  is  due  to  the  presence  of  seawater  scat¬ 
tering  for  r  =  0.05  m  and  p  =  0.99. 


The  analytical  forms  of  the  sensitivity  coefficients 
derived  in  this  section  are  valid  for  any  PSICAM 
application;  however,  the  ranges  of  the  physical  pa¬ 
rameters  in  our  examples,  figures,  and  discussions 
are  chosen  to  be  consistent  with  applications  to 
coastal  ocean  waters,  visible  wavelengths,  and  a  PSI¬ 
CAM  cavity  made  with  Spectralon  (Labsphere,  Inc., 
North  Sutton,  N.H.).  The  value  of  the  absorption 
coefficient  of  pure  seawater  for  light  in  the  400- 
700-nm  range  has  a  minimum  value  of  0.014  m  1  (at 
430  nm)  and  a  maximum  value  of  0.65  nW1  (at  700 
nm).  The  value  of  a  in  this  wavelength  range  for 
natural  ocean  waters  is  typically  between  0.1  and  1.0 
m  ',  depending  on  the  wavelength  and  the  concen¬ 
trations  of  phytoplankton,  suspended  sediments,  and 
suspended  and  dissolved  organic  material.  The  re¬ 
flectivity  of  Spectralon  is  approximately  p  =  0.99. 

A.  Known  p 

If  the  values  of  p,  r,  and  aref  are  known,  we  can 
determine  the  absorption  coefficient  of  a  fluid  sample 
by  measuring  T  [Eq.  (1)]  and  then  numerically  solv¬ 
ing  Eq.  (10).  Here  we  examine  the  sensitivity  of  the 
value  of  a  obtained  from  Eq.  (10)  to  the  four  quanti¬ 
ties  T,  p,  r,  and  aref.  In  general  the  values  of  the  four 
sensitivity  coefficients  Ba/BT,  Ba/Bp,  3a/ Br,  and  Ba/ 
Bare{  depend  on  the  four  physical  parameters  p,  r,  aref, 
and  a. 

For  T  given  by  Eq.  (10)  and  Ps  given  by  Eq.  (6), 

T  Ba  _ 
a  BT 

2a2r2  -  p[l  -  (2 ar  +  l)exp(-2ar)] 

2a3r3  -  p[(2a2r2  +  3ar  +  2)exp(-2ar)  -  2  +  ar]  ’ 

(11) 

where  we  have  ignored  the  size  of  the  source.  Note 
that  on  the  right-hand  side  of  Eq.  (11)  the  parameters 
a  and  r  exist  only  as  the  product  (ar).  We  show  the 
normalized  form  of  the  sensitivity  coefficient  in  Eq. 
(11)  because  it  is  conveniently  independent  of  aref. 
We  can  find  the  value  of  Ba/ BT  by  multiplying  Eq.  (11) 
by  (a/71)  with  T  replaced  by  the  expression  in  Eq.  (10), 
which  introduces  a  dependence  on  aref.  The  sensi¬ 
tivity  coefficient  of  a  with  respect  to  p  can  be  found 
from  Eq.  (10)  with 


Ba 

BT/Bp 

Bp 

BT/Ba ' 

Likewise 

Ba  BT/Br 

Ba 

BT/Ba„  f 

Br  BT/Ba  ’ 

Bare  f 

BT/Ba 

Table  1  gives  the  values  of  the  four  sensitivity  co¬ 
efficients  in  both  their  unnormalized  and  their  nor¬ 
malized  form  for  the  example  cases  of  a  =  0.2  m_1, 
a  ref  =  0.02  m_1,  r  =  0.05  m,  and  p  =  0.99.  We 
verified  these  coefficients  with  maple  v  (Waterloo  Ma¬ 
ple,  Inc.,  Waterloo,  Ontario)  by  solving  Eq.  (10)  with 
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Table  1.  Sensitivity  Coefficients  for  the  Determination  of  a  with 
Eqs.  (6)  and  (10)  when  p  is  Known" 


Coefficient 

Value 

da/ dT 

-0.72  nV1 

da /dp 

-16  m-1 

da /dr 

—3.2  m~2 

da/dare  f 

2.1 

(T/a)(da/dT) 

-1.7 

(p/a)(da/dp) 

-79 

(r/a)(da/dr) 

-0.79 

(®ref/ a)(da/ Balet) 

0.21 

“ a  =  0.2  m  1,  aref  =  0.02  m  1,  r  =  0.05  m,  and  p  =  0.99. 


and  without  a  small  error  in  each  parameter.  Each 
unnormalized  sensitivity  coefficient  gives  the  ratio  of 
the  absolute  error  in  a  to  the  absolute  error  in  the 
parameter  causing  the  error  in  a,  whereas  each  nor¬ 
malized  sensitivity  coefficient  gives  the  ratio  of  the 
percentage  error  in  a  to  the  percentage  error  in  the 
given  parameter.  It  can  be  seen  from  Table  1  that 
the  value  of  a  is  most  sensitive  to  p;  for  the  parameter 
values  used  in  Table  1,  a  1%  error  in  p  would  cause 
nearly  an  80%  error  in  a.  The  values  of  the  sensi¬ 
tivity  coefficients  of  a  with  respect  to  T,  aref,  and  r  are 
much  smaller  but  nonetheless  potentially  significant. 
Letting  AT  represent  the  uncertainty  in  the  measure¬ 
ment  of  T  and  Ap  represent  the  uncertainty  in  the 
measurement  of  p,  etc.,  the  total  uncertainty  in  our 
computed  value  of  a  is12 


Aa 


da 


+ 


dT 
da 

d^ref, 


da 

dp 

1/2 


(AT)  +  —  (Apr  + 


(Ar)2 


(Aaref)2 


(14) 


If,  for  the  example  used  in  Table  1,  we  measure  r 
(0.050  m)  to  ±0.1  mm,  p  (0.990)  to  ±0.002,  aref  (0.020 
m_1)  to  ±0.005  m~  ,  and  the  corresponding  value  of 
T  (0.486)  to  1%  (i.e.,  F0  is  measured  to  ±0.5%),  then 
from  Eq.  (14)  and  Table  1 


A  a  =  {[(0.71)(0.00486)]2  +  [(16)(0.002)]2  +  [(3.2) 

X  (0.0001)]2  +  [(2.1)(0.005)]2}1/2  m  1 
=  (1.2  X  1CT5  +  1.0  X  1(T3  +  1.0  X  1(T7 

+  1.1  X  10-4)1/2  nT1 

=  0.034  nT1.  (15) 


Therefore  our  computed  value  of  a  would  be  0.200  ± 
0.049  m-1,  with  the  uncertainty  due  primarily  to  the 
measurements  of  p  and  aref  and  the  uncertainty  in  r 
being  insignificant. 

We  are  perhaps  more  interested  in  the  absolute 
uncertainty  in  a  that  is  due  to  the  relative  (percent¬ 
age)  uncertainty  in  T,  which  is  quantified  by  the  term 
(T)(da/dT),  than  in  either  ( da/dT )  or  (T/a)(da/dT). 
We  can  find  iT)( da/dT)  by  simply  multiplying  Eq.  (11) 
by  a,  and  the  term  is  shown  versus  a  in  Fig.  4.  It  is 
negative,  indicating  that  the  value  of  a  will  be  under- 


Fig.  4.  Seminormalized  sensitivity  coefficient  (T)(da/dT)  (m  1) 
versus  a  (m_1). 


estimated  if  the  measured  value  of  T  is  greater  than 
the  true  value  and  vice  versa.  It  can  be  seen  from 
Fig.  4  that  |T(3a/dT)|  increases  approximately  lin¬ 
early  with  increasing  values  of  a  and  decreases  with 
increasing  values  of  r  and  p.  It  was  found  that 
\T{da/dT}\  is  relatively  insensitive  to  the  value  of  r  for 
r  greater  than  approximately  0.02  but  increases  rap¬ 
idly  with  decreasing  r  for  r  <  0.02.  The  increase  in 
\T(da/dT)\  with  decreasing  p  is  approximately  linear. 

The  value  of  da/ dp  (Fig.  5)  equals  zero  when  a  =  aref 
and  decreases  approximately  linearly  with  increas¬ 
ing  a.  Increasing  aref  or  decreasing  p  moves  the 
slope  of  the  (da /dp)  curve  versus  the  a  curve  toward 
zero.  Increasing  r  also  moves  the  slope  toward  zero; 
however,  the  effect  is  small  and  is  not  illustrated  in 
Fig.  5.  Similarly  the  value  of  da/ dr  (Fig.  6)  equals 
zero  when  a  =  aref  and  decreases  approximately  lin¬ 
early  with  increasing  a.  The  slope  of  (da/ dr)  versus 
a  is  steepest  for  small  values  of  aref.  The  value  of 
da/daTef  (Fig.  7)  equals  unity  when  a  =  aref,  increases 
approximately  linearly  with  a,  and  is  largest  for 
small  values  of  aref.  It  was  found  that  the  slope  of 
the  (da/daref)  curve  versus  the  a  curve  increases  when 
the  value  of  p  or  r  increases. 

It  is  important  to  note  from  Figs.  5-7  that  the  value 
of  arefhas  a  dramatic  effect  on  the  overall  uncertainty 


a 

Fig.  5.  Sensitivity  coefficient  da/d p  (m_1)  versus  a  (m_1)  for  r  = 
0.05  m. 
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Fig.  6.  Sensitivity  coefficient  da/dr  (m  2)  versus  a  (m  4)  when  p 
is  known  for  r  =  0.05  m  and  p  =  0.99. 

in  a.  In  particular,  we  can  reduce  both  sensitivity 
coefficients  da/d p  (Fig.  5)  and  da /dr  (Fig.  6)  to  values 
near  zero  by  taking  aref  *=  a.  Although  the  values  of 
r  and  p  are  generally  fixed  once  a  PSICAM  is  con¬ 
structed,  we  can  control  the  value  of  aref  by  varying 
the  concentration  of  a  nonscattering  dye.  We  could 
obtain  the  most  precise  determination  of  a  by  first 
measuring  F0  for  the  unknown  sample  and  then  by 
adjusting  the  dye  concentration  in  the  reference  so¬ 
lution  as  many  times  as  necessary  until  T  =  1  (i.e., 
aref  =  a).  In  this  way  a  reference  solution  would  be 
produced  that  has  the  same  absorption  coefficient  as 
the  sample  of  interest  but  presumably  has  a  much 
smaller  scattering  coefficient.  We  could  then  mea¬ 
sure  the  value  of  a  =  aref  by  placing  the  reference 
solution  in  a  spectrophotometer.  The  total  uncer¬ 
tainty  in  a  would  be  approximately  equal  to  that  in 
aref  (3a/daref  =  1  when  aref  =  a)  plus  that  due  to  the 
uncertainty  in  the  measurement  of  T[T(da  /  dT)  is  in¬ 
dependent  of  the  value  of  aref].  If  aref  in  the  example 
in  Eq.  (15)  were  changed  to  aref  =  a,  the  uncertainty 
in  a  would  be  reduced  from  0.034  m  1  to 

A  a  =  {[(0.71)(0.00486)]2  +  0  +  0  +  [(1) 

X  (0.005)]2}172  m  1  =  0.0061  m"1.  (16) 


Fig.  7.  Sensitivity  coefficient  da/daTef  when  p  is  known  for  r  = 
0.05  m  and  p  =  0.99. 


Table  2.  Sensitivity  Coefficients  for  the  Determination  of  a  with  Eqs. 
(6)— (10)  when  p  is  Unknown" 


Coefficient 

Value 

da/dT 

-0.72  m^1 

da/Baref  ( B  +  R ) 

2.1 

da/ dTjyz 

0.83  m^1 

da/daA 

4.2 

da/daB  ( B  =£  R) 

-5.3 

da/dr 

-0.0026  m~2 

da/daref  (B  =  R) 

-1.2 

C DOa/dT ) 

-0.35  m-1 

(alef/a)(da/daret) 

0.21 

(Tab)^0/ d Tab ) 

0.56  m-1 

(aA/a)(da/daA) 

2.1 

(aB/a)(da/daB) 

-2.6 

(. r/a)(da/dr •) 

-6.5  X  10-4 

(are{/a)(da/daTJ 

-0.12 

“a  =  0.2  m  1,aref=  0.02  m  1,  r  =  0.05  m,  p  =  0.99,  aA  =  0.1  m  1, 
and  aB  =  0.05  m-1. 


Admittedly,  however,  it  may  be  prohibitively  difficult 
to  produce  many  reference  solutions.  Also  note  that, 
although  da/d p  is  theoretically  zero  when  a  =  aref,  the 
slope  of  da/d p  versus  a  (Fig.  5)  is  in  some  cases  so 
steep  that  the  uncertainty  in  aref  introduces  a  signif¬ 
icant  uncertainty  in  da/d p. 

B.  Unknown  p 

If  the  value  of  p  to  be  used  in  Eq.  (10)  is  unknown  it 
can  be  determined  with  Eq.  (9).  In  this  case  our 
estimate  of  a  ultimately  depends  on  the  measure¬ 
ments  of  r,  aA,  aB,  TAB,  aref,  and  T.  The  value  of  the 
sensitivity  coefficient  da/dT  is  the  same  as  that  given 
in  Subsection  4.A  [Eq.  (11)  and  Fig.  4].  The  value  of 
da/daTe{,  on  the  other  hand,  is  the  same  as  that  given 
in  Subsection  4.A  only  if  the  reference  sample  B  used 
in  the  determination  of  p  [Eq.  (9)]  is  not  chosen  to  be 
the  same  as  the  reference  sample  R  used  in  the  de¬ 
termination  of  a  [Eq.  (10)].  Because  the  measure¬ 
ment  Tab  depends  on  r,  the  value  of  da/ dr  is  not  the 
same  as  that  given  in  Subsection  4.A.  We  therefore 
need  to  determine  the  new  sensitivity  coefficients  da/ 
dr,  da/daA,  da/daB,  da/dT AR,  and,  when  R  and  B  are 
the  same  fluid,  da/dan.f.  Each  new  sensitivity  coef¬ 
ficient  depends  in  general  on  the  six  physical  param¬ 
eters  r,  p,  aA,  aB,  aref,  and  a. 

We  know  from  Subsection  4.A  that  for  the  special 
case  of  aref  =  a  the  values  of  p  and  r  in  Eq.  (10)  have 
no  effect  on  the  estimate  of  a.  Because  the  values  of 
aA,  aB,  and  Tab  affect  a  only  through  p,  these  three 
measurements  are  unimportant,  and  therefore  inef¬ 
fective,  when  aref  =  a.  Therefore  we  no  longer  want 
to  set  aref  a  when  using  Eqs.  (9)  and  (10)  in  com¬ 
bination.  Instead  we  take  advantage  of  the  fact  that 
Eq.  (9)  is  the  inverse  of  Eq.  (10).  The  measurement 
error  in  r  tends  to  cancel  when  these  equations  are 
used  in  combination.  Likewise,  if  we  take  R  and  B 
to  be  the  same  solution,  then  the  estimate  of  a  will  be 
relatively  insensitive  to  the  measurement  of  aref. 

Table  2  gives  the  values  of  the  sensitivity  coeffi¬ 
cients,  both  in  their  unnormalized  and  in  their  nor- 
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malized  form,  for  the  case  of  a  =  0.2  m-1,  aref  =  0.02 
nT1,  aA  =  0.1  nT1,  aB  =  0.05  nT1,  r  =  0.05  m,  and 
p  =  0.99.  Note  that  all  the  sensitivity  coefficients 
have  relatively  small  values  compared  with  that  of 
da/d p  from  Eq.  (10)  alone  (Table  1).  Also  note  in 
Table  2  that  da/dr  is  now  -0.0026  m  2,  compared 
with  — 3.2  nT2  in  Table  1,  verifying  that  the  precision 
of  the  measurement  of  r  is  even  less  important  when 
a  is  determined  with  Eqs.  (9)  and  (10)  together  than 
when  a  is  found  from  Eq.  (10)  alone.  (The  uncer¬ 
tainty  of  r  is  unlikely  to  be  important  in  either  case.) 
Furthermore,  we  can  see  from  Table  2  that  to/toref 
when  B  =  R  is  smaller  than  either  da/ d aref  or  da/ daB 
when  B  A  R,  confirming  that  we  should  let  B  and  R 
be  the  same  solution.  Accordingly,  we  assume  that 
B  =  R  for  the  remainder  of  this  subsection.  Ignoring 
Ar,  we  obtain  the  uncertainty  in  a  from 


Aa  = 


,  to  Wat 

dT 


+  T 


to 


dT, 


AB, 


+ 


AT 


to 

A^ref, 
2 


(Aaref)2 


AB 


[  AB 


+W(i->2 


1/2 


(17) 


If,  for  the  example  used  in  Table  2,  we  measure  aref 
(0.02  m-1)  and  aA  (0.1  m-1)  each  to  ±0.005  nT1  and 
the  corresponding  values  of  T  and  Wb  each  to  ±1%, 
then  from  Eq.  (17)  and  Table  2  we  obtain 

A  a  =  {[(0.35)(0.01)]2  +  [(1.2)(0.005)]2  +  [(0.56) 

X  (0.01)]2  +  [(4.2)(0.005)]2}1/2  nT1 
=  (1.2  X  10~5  +  3.6  X  1CT5  +  3.1  X  10“5 
+  4.4  X  10W1/2  nT1  =  0.023  nT1.  (18) 

Here  the  uncertainty  in  a  is  due  mostly  to  the  uncer¬ 
tainty  in  aA.  The  calculated  value  of  a  (0.200  ± 
0.036  nT1)  is  an  improvement  over  the  example 
given  in  Eq.  (15),  indicating  that  even  if  the  value  of 
p  is  known  from  direct  measurement  it  may  be  pru¬ 
dent  to  use  Eq.  (9)  to  determine  the  value  of  p  for  use 
in  Eq.  (10)  rather  than  the  measured  value  of  p.  Be¬ 
cause  the  experimental  setup  is  the  same  for  Eqs.  (9) 
and  (10),  the  uncertainties  in  the  physical  parame¬ 
ters  are  far  less  important  when  a  combination  of 
Eqs.  (9)  and  (10)  is  used  than  when  Eq.  (10)  is  used 
alone.  Even  if  the  directly  measured  value  of  p  is 
more  accurate  than  the  value  of  p  obtained  from  Eq. 
(9),  it  generally  will  not  lead  to  a  more  accurate  value 
of  a  when  used  in  Eq.  (10). 

The  derivation  of  general  expressions  for  the  sen¬ 
sitivity  coefficients  for  Eqs.  (9)  and  (10)  is  outlined  in 
Appendix  B.  The  value  of  to /dr  was  found  to  be 
negligibly  small  over  a  wide  range  of  physical  param¬ 
eters  and  is  identically  zero  when  a  =  aA.  Although 
an  error  in  r  results  in  an  error  in  p  with  Eq.  (9),  the 
errors  in  r  and  p  are  such  that  their  effects  cancel 
each  other  in  the  determination  of  a  with  Eq.  (10). 
The  value  of  to/toref  is  also  zero  when  aA  =  a. 
Therefore  the  estimate  of  a  from  Eqs.  (9)  and  (10) 


Fig.  8.  Sensitivity  coefficient  da/daA  for  r  =  0.05  m  and  p  =  0.99. 


depends  only  on  the  measurements  of  T,  TAB,  and  aA 
when  aA  «  a. 

As  illustrated  in  Fig.  8  the  value  of  da/daA  is  ap¬ 
proximately  linearly  dependent  on  a  and  equals  zero 
when  a  =  aref.  The  slope  of  the  da/ da A  curve  versus 
the  a  curve  is  larger  for  small  values  of  aA  (i.e.,  aA  > 
a)  than  for  large  values  of  aA.  It  is  positive  when 
aA  >  aref  and  negative  when  aA  <  aref.  Similarly,  as 
shown  in  Fig.  9,  the  value  of  TAB(da/dTAB)  is  approx¬ 
imately  linearly  dependent  on  a  and  equals  zero 
when  a  =  aref.  The  slope  of  (TABda/dTA B)  versus  a  is 
positive  when  aA  >  aref  and  negative  when  aA  <  aref. 

We  conclude  that  it  is  desirable  for  aA  to  be  at  least 
as  large  as  a.  For  aA  >  a,  we  are  faced  with  a  trade-off: 
increasing  aA  generally  decreases  the  sensitivity  of  a 
to  both  aA  and  TAB  but  increases  its  sensitivity  to  aref 

5.  Nondiffuse  Source 

Up  to  this  point  in  our  analysis  we  have  assumed  that 
the  PSICAM  is  illuminated  with  a  perfectly  isotropic 
source.  However,  a  small  isotropic  source  is  techni¬ 
cally  difficult  to  construct.  Therefore  in  this  section 
we  evaluate  the  performance  of  a  PSICAM  when  the 
source  is  anisotropic.  Three-dimensional  Monte 


Fig.  9.  Seminormalized  sensitivity  coefficient  (T^ida / dTj^)  ver¬ 
sus  a  (m_1)  for  r  =  0.05  m  and  p  =  0.99. 
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Fig.  10.  Percentage  error  in  the  PSICAM  detector  response  as  a 
function  of  detector  location  for  anisotropic  light  sources.  The 
detector  position  is  given  by  angle  0rf  with  respect  to  the  forward 
direction  of  the  source.  The  source  angular  distribution  is  the 
Henyey-Greenstein  function  with  asymmetry  factors  g  =  0,  0.2, 
and  0.5.  The  results  shown  are  for  a  =  0.3  m_1,  b  =  0.7  m_1,  r  = 
0.05  m,  and  p  =  0.99. 


Carlo  simulations  were  performed  for  several  differ¬ 
ent  types  of  centrally  located  source,  and  the  wall 
irradiance  was  determined  at  various  sections  of  the 
PSICAM  wall.  It  is  important  to  note  that  the  total 
irradiance  incident  on  the  entire  cavity  wall  is  inde¬ 
pendent  of  the  shape  of  the  source  function;  however, 
the  irradiance  at  a  particular  location  on  the  wall, 
and  therefore  at  the  detector,  does  depend  on  the 
source  shape. 

Shown  in  Fig.  10  is  the  percentage  error  in  PSI¬ 
CAM  detector  irradiance  for  central  light  sources 
that  are  characterized  by  the  Henyey-Greenstein 
function  [analogous  to  Eq.  (A4)]  for  asymmetry  factor 
values  oig  =  0,  0.2,  and  0.5  and  for  r  =  0.05  m,  p  = 
0.99,  a  =  0.3  m-1,  and  b  =  0.7  m-1.  The  error  is 
plotted  against  the  position  of  the  detector,  expressed 
in  terms  of  the  angle  0d  between  the  forward  direction 
of  the  light  source  and  the  radius  from  the  center  to 
the  detector.  (For  example,  0d  =  0  when  the  source 
is  pointed  straight  at  the  detector.)  If  the  source  is 
perfectly  isotropic,  the  irradiance  is  identical  at  all 
positions  on  the  cavity  wall  and  the  error  that  is  due 
to  the  source  is  zero  everywhere.  For  a  source  char¬ 
acterized  by  the  Henyey-Greenstein  function,  the  ir¬ 
radiance  is  stronger  than  that  for  an  isotropic  source 
at  small  values  of  0d,  whereas  the  irradiance  is 
weaker  than  that  for  an  isotropic  source  at  large  val¬ 
ues  of  0d.  For  a  source  with  g  =  0.2,  65%  of  the  total 
source  light  is  projected  in  the  forward  180  deg,  most 
of  which  is  concentrated  in  the  near-forward  direc¬ 
tions,  and  35%  is  projected  in  the  backward  180  deg. 
However,  because  of  the  large  number  of  wall  reflec¬ 
tions,  in  this  case  33,  the  error  in  the  wall  irradiance 
is  less  than  2.5%  at  all  locations.  It  can  be  seen  from 
Fig.  10  that  it  would  be  best  to  place  the  detector  at 
approximately  a  right  angle  to  the  sensor,  where  the 
error  that  is  due  to  the  source  is  near  zero.  For  a 
source  with  g  =  0.5,  83%  of  the  source  light  is 


projected  in  the  forward  180  deg  and  17%  in  the 
backward  directions.  Although  this  extremely 
anisotropic  source  would  not  be  practical,  we  con¬ 
clude  from  Fig.  10  that,  theoretically,  it  could  be  used 
if  the  detector  were  placed  near  0d  =  60°. 

It  is  important  to  note  that  the  error  in  T  or  Tab 
may  be  much  smaller  than  the  error  in  the  detector 
response.  If  the  detector  receives  more  source  light 
than  it  would  with  a  perfectly  isotropic  source,  it  will 
do  so  both  with  the  reference  fluid  and  with  the  sam¬ 
ple  of  interest.  The  error  will  not  completely  cancel 
out,  however,  because  the  amount  of  scattering  and 
the  average  number  of  wall  reflections  will  be  differ¬ 
ent  in  the  reference  than  in  the  sample. 

In  summary  we  conclude  that,  although  it  is  not 
critical  to  construct  the  light  source  to  be  perfectly 
isotropic,  the  error  that  is  due  to  this  anisotropy 
should  be  minimized  when  the  irradiance  sensor  is 
positioned  at  a  location  where  the  source  strength  is 
near  its  average  value. 

If  a  sufficiently  isotropic  source  cannot  be  obtained, 
it  may  be  preferable  to  replace  the  plane  irradiance 
sensor  with  a  radiance  sensor  that  is  positioned  such 
that  the  source  is  not  within  the  detector’s  field  of 
view.  In  this  case  Eq.  (2)  would  be  replaced  with 


Nc  = 


pPpPs 

1  -  pPs  ’ 


(19) 


an  approximation  that  ignores  the  scattering  of  light 
during  its  first  pass  from  the  source  to  the  wall.  We 
will  quantify  the  accuracy  of  Eq.  (19)  in  a  future  paper. 


6.  Conclusions 

Consistent  with  Kirk,9-10  we  found  that  the  scattering 
errors  in  Eqs.  (4)  and  (6)  are  insignificant  for  most 
ocean  optics  problems.  However,  as  demonstrated 
in  Figs.  1-3,  the  scattering  errors  may  be  important 
for  near-IR  light  or  for  other  PSICAM  applications. 

If  the  value  of  the  PSICAM  wall  reflectivity  is 
known,  the  absorption  coefficient  of  a  sample  can  be 
determined  [with  Eq.  (10)]  from  PSICAM  measure¬ 
ments  of  the  sample  and  of  a  reference  solution.  Be¬ 
cause  of  the  high  sensitivity  of  a  to  p,  however,  the 
uncertainty  in  the  result  may  be  unacceptably  large 
unless  p  is  known  to  high  precision  and  a  aref.  On 
the  other  hand,  we  can  obtain  a  precise  determina¬ 
tion  of  a  by  producing  a  nonscattering  dye  solution 
such  that  T  =  1  and  by  determining  the  value  of  aref 
with  a  spectrophotometer. 

Alternatively,  we  can  obtain  a  generally  more  pre¬ 
cise  result  by  including  the  PSICAM  measurement  of 
at  least  one  additional  reference  solution,  eliminating 
the  need  to  know  the  value  of  the  wall  reflectivity  a 
priori.  Even  if  the  value  of  p  is  known  from  direct 
measurement,  it  may  be  preferable  to  determine  it 
with  Eq.  (9)  because  this  tends  to  eliminate  errors 
associated  with  the  measurements  of  p,  r,  and  aref. 

A  prototype  PSICAM  is  currently  under  construc¬ 
tion,  and  we  intend  to  publish  a  description  of  the 
instrument  design  and  laboratory  experiments  once 
they  have  been  completed.  We  will  also  present  an 
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analysis  of  a  PSICAM  in  which  a  radiance  sensor  is 
used  [with  Eq.  (19)]  in  place  of  a  plane  irradiance 
sensor.  A  radiance  sensor  is  easier  to  implement 
and  may  reduce  errors  associated  with  a  nonisotropic 
light  source;  however,  it  also  introduces  a  new  source 
of  scattering  error.  The  significance  of  the  scatter¬ 
ing  error  must  be  evaluated  for  the  range  of  ocean 
water  types  before  one  can  safely  adopt  the  radiance 
approach. 

Appendix  A:  Monte  Carlo  Techniques 

In  the  Monte  Carlo  simulations  photons  were  traced 
one  at  a  time  and  the  light  field  was  taken  to  be  the 
superposition  of  all  the  photon  paths.  The  path  of 
each  photon  was  determined  by  random  numbers  91 
evenly  distributed  in  the  range  [0,  1].  The  photon 
path  length  s  was  found  from 


sines  of  the  photon  (a',  (3',  y')  were  computed  from  the 
initial  photon  direction  cosines  (ct,  (3,  y)  and  scatter¬ 
ing  angles  (polar  angle  ©,  azimuthal  angle  <h )  with 


a' 

p' 

— 

ay/  Vl  -  y2 
|3y/  Vl  -  y2 

-(3/  Vi  -  y2  a 
a/  Vl  -  y2  P 

Y_ 

<N 

1 

rH 

1 

_ 1 

0  y 

X 


,  (y2  <  1), 


a' 

as 

P' 

=  sign(y) 

ps 

j'. 

_ys_ 

,  (72  =  D, 


where 


(A8) 


(A9) 


s  =  — (l/c)ln(l  -  91),  (Al) 

where  c  is  the  attenuation  coefficient  (1/m).  At  each 
scattering  event,  the  azimuthal  scattering  angle  with 
respect  to  the  incident  direction  was  found  with 

$  =  2ir9t.  (A2) 

For  isotropic  scattering,  the  cosine  of  the  polar  scat¬ 
tering  angle  with  respect  to  the  incident  direction  |jls. 
was  found  with 

(xs  =  1  -  291.  (A3) 

For  Henyey-Greenstein  scattering13-14 


as  =  sin(@)cos(<l>),  (3S  =  sin(©)sin(<b),  ys  =  cos(©). 

(A10) 

Appendix  B:  Sensitivity  Coefficients  when  p  is  Known 

Let  T4  designate  T{a,  aref,  r,  p)  as  given  by  Eq.  (10): 

T4  =  T(a,  aref,  r ,  p).  (Bl) 

The  derivatives  da/dT  and  da/daref  (B  #  R)  are  the 
same  as  those  given  in  Subsection  4.A: 

da  1 1  c)T4\ 

—  <a,«W,p)  =  l/]— ),  (B2) 


P(g;  pJ 


i  i-g2 
4tt  (1  +  g2  ~  2g-p.s)3/2  ’ 


(A4) 


da 

d&ref 


(a,  aref,  r,  p) 


,  (B  +  R). 


(B3) 


where  g  is  the  scattering  asymmetry  factor.  The 
corresponding  cumulative  distribution  function  is 


C(ps)  =  ; 


i  -g2 


2  I  (1  +  g2  -  2g^s)3/2 


i-g2 

2g 


l~g  (1  +g2  -  2£|xs)1/2 


(A5) 


Setting  91  equal  to  the  cumulative  distribution  func¬ 
tion  and  solving  for  |xs,  we  found  jjis.  with 


The  other  sensitivity  coefficients,  however,  depend  on 
six  parameters:  interchangeably  (a,  aref,  r,  aA,  aB, 
Tab )  or  (a,  aref,  r,  aA,  aB,  p).  Let  p(aA,  aB,  r,  T^)  be 
given  by  Eq.  (9)  and  aB,  r,  p)  be  given  by  Eq. 

(3).  We  can  substitute  p(aA,  aB,  r,  TAB)  into  the 
expression  for  T(a,  aref,  r,  p)  and  thereby  express  T  in 
terms  of  six  parameters: 

Tq  =  T(a,  aie f,  r ,  aA,  aB,  TA/1).  (B4) 

Each  of  the  sensitivity  coefficients  can  be  determined 


2g  +  1  -  291g-  -  291  +  2912g  +  g2  -  2g3(3l  -  2°J\g2  +  2§3912 
i-g~  1  +  291 gf 


(A6) 


For  Petzold  scattering,  we  solved  for  p,s  with 


C(ps)  =  91,  (A7) 

where  the  cumulative  distribution  function  Cfp.,)  was 
taken  to  be  the  average  of  those  provided  by  Petzold11 
for  the  three  samples  from  San  Diego  Harbor. 

For  each  scattering  event,  the  new  direction  co- 


either  directly  from  T6  or  from  a  combination  of  T A 
and  p  as  follows: 


dT  dT6 

—  (a,  aref,  r,  aA,  aB,  p)  =  — 
dr  dr 

dT4  dT4  3p 

= - + - -, 

dr  3p  dr 


(B5) 
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—  (a,  a  reft  r,  aA,  aB,  p)  =  -(dTJ  dr) /(dTJ  da ) 
dr 


c)T4  dp 
dp  dr 


(dTJ  (da), 


(B6) 


—  (a,  a  reft  r,  aA,  aB,  p)  =  -( dTJ  daA)  /  (dTJ  da ) 
daA 

=  -[(dT4/dp)(dp/daA)]/(dT4/da),  (B7) 


—  (a,  aref,  r,  aA,  aB,  p)  =  ~(dTJdaB)/(dTJda) 
daB 

=  ~\{dTJ  dp)(dp/ daB)y  {dTJ  da),  (B8) 


(a,  aref,  r,  aA,  aB,  p)  =  -(dTJ dTAB)/ (dTJ da) 

dl-AB 

=  -[(dTJdp)(dp/dTAB)]/(dTJda),  (B9) 

If  the  same  reference  solution  is  used  in  Eqs.  (9) 
and  (10)  (i.e.,  B  =  R ),  then  T  can  be  expressed  in 
terms  of  five  parameters: 


T5  —  T(a,  a  ref,  r ,  ctA,  Tab),  (BIO) 

and  the  sensitivity  coefficient  for  aref  becomes 


da 

daref 


(a,  aref,  r,  aA,  p) 


9T 4  |  /ar4V  dp  \ 
daref  \dp)\daiefJ 


(Bll) 
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